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1.  Introduction 

The  computation  of  probabilities  associated  with  sequential  stopping  variables 
can  often  be  easily  performed  by  recursive  techniques  (e.g.  Armitage  et  al  ( 1969), 
Aroian  and  Robinson  (1969),  Samuel-Cahn  (1974),  Pocock  (1977),  Zacks  (1974,  1980, 
1981,  1987)).  Explicit  formulae  for  probabilities  of  boundary  crossings  are  generally 
difficult  to  obtain.  Wiener  process  and  other  asymptotic  approximations  to  the 
distributions  of  stopping  times  are  discussed  by  Siegmund  (1985).  The  present  paper 
shows  that  in  the  case  of  sequential  procedures  based  on  Poisson  processes,  explicit 
expression  can  be  derived.  These  derivations  are  facilitated  by  the  nondecreasing 
nature  of  the  sample  paths,  the  memory-less  property  of  the  Poisson  process,  and 
the  possibility  to  solve  certain  recursive  relations  associated  with  the  probability 
distributions  of  the  stopping  variables.  For  the  purpose  of  showing  the  potential  of 
this  technique,  we  derive  explicit  expressions  for  the  cumulative  distribution 
functions  (CDF)  of  stopping  variables  associated  with  the  acceptance  and  with  the 
rejection  boundaries  of  a  Wald’s  SPRT,  which  is  frequency  truncated.  Such 
stopping  boundaries  are  applied  in  reliability  life  testing  for  systems  having 
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exponentially  distributed  life  length  (see  MIL-STD  781C  (1977)).  Methods  of 
determining  confidence  intervals  for  the  MTBF,  after  such  sequential  life  testing, 
were  given  by  Siegmund  (1978),  Bryant  and  Schmee  (1979),  and  Zacks  (1987). 

Siegmund  approximated  the  distributions  of  stopping  times,  while  Zacks  determined 
them  recursively.  The  results  of  the  present  paper  eliminate  the  need  for 
recursive  computations  in  the  Poisson  case. 

In  Section  2  we  present  the  stopping  boundaries,  the  observed  process  and 
some  general  definitions.  In  Section  3  we  deriva  explicitly  certain  boundary  crossing 
probabilities.  Section  4  is  davoted  to  an  actual  derivation  of  the  distributions  of 
the  stopping  variables,  for  a  particular  case  under  consideration.  A  numerical 
example  is  provided  in  Section  5. 


2.  The  Failure  Process  and  Stopping  Variables 

Let  {X(tk  9  <  t  <  oo)  be  a  Poisson  process  with  intensity  X,  0  <  X  <  o©;  i.eM 
E^(X(t)>  —  Xt,  0  <  t  <  oo.  We  consider  a  truncated  Wald  SPRT,  which  is  specified 
by  the  boundary  lines 

bt,  0  <  t  £  tkr.kj 

Vk^  1  *  Vk() 

and 


(2.1) 


bR(t) 


-( 


k,+ 


(2.2) 


bA(t) 


-  kn  4-  bt, 


£  t  £  t 


where  kj,  kQ,  ks  are  integers,  ks  >  kj,  and  tj  «=  i/b  (1  «=  1,  2,  ).  The  boundary 
b^(t)  is  called  the  acceptance  boundary,  and  bR(t)  is  called  the  rejection  boundary. 
Notice  that  b^(t)  and  bR(t)  intersect  at  (t^  +k  ,  ks). 

0  5 

We  define  the  stopping  times  (variables) 


(2.3) 


tR 

ta 


inf  {t:  X(t)  ^  bR(t),  0  <  t  <  \ 

inf {t:  X(t)  =  bA(t),  tkQ  <;  t  <.  tko+kg}, 


(2.4) 
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where  |X(t):  t  >  oj  is  a  nondecreasing  jump  process,  with  jumps  of  one  unit.  The 
process  is  stopped  at  the  first  instant  r  at  which  it  crosses  either  the  lower 
boundary,  bA(t),  or  the  upper  boundary  bjj(t).  Obviously,  t  =■=  min  <rA,  r^h 

For  the  derivation  of  the  distributions  of  the  stopping  times  t^,  ta  and  r 
one  has  to  distinguish  between  three  cases, 

Case  I:  kj<  ks  £  kj  4-  kQ; 

Case  II:  kj  +  kg<  ks  £  2kj  +  kQ; 

Case  III:  2k j  -f  kg<  ks. 

Since  the  method  of  derivation  of  the  distributions  is  similar  in  all  the  three  cases, 
we  focus  attention  in  the  present  paper  only  on  Case  I. 

The  sample  paths  of  the  Poisson  process  {X(t>i  0  <  t)  assume  only 
nonnegative  integer  values.  Thus,  a  sample  path  can  cross  the  lower  boundary  bA(t) 
only  at  the  discrete  times  tj,  kQ  £  i  ^  k8  +  kg  —  1.  Let 


C0 « - 


tk0+j) 


be  the  probability  distribution  function  of  the  discrete  stopping  variable  rA,  over 
tkg+j.  We  introduce  the  notation 

{1,  if  A  is  true 

0,  otherwise; 

and  p(j;0  and  Pos(j;£)  are  to  be,  respectively,  the  PDF  and  CDF  of  the  Poisson 
distribution  with  mean  $. 

(X) 

Let  F^  (t)  denote  the  corresponding  CDF  of  r^.  The  stopping  time  on 
the  other  hand,  has  an  absolutely  continuous  distribution.  Define  the  probability 
f  unction 


gx(j;  t)  =  Px{X(t)  =  j,  t  ^  t},  j  =  0,  1, ... 


WWW 


mj  I. 
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Let  [a]  be  the  largest  integer  smaller  than  or  equal  to  a,  and  define 
,  (kj  4  btj,  if  kj  •+-  bt  is  not  an  integer 


kj4-  bt  -  1, 


otherwise. 


Similarly,  let  a"1*  «=  max  {a,  0}  and  let 

[-  kQ  4-  bt], 

(2.7)  m(t)  —  I 

l  [-  k 


if  t  ^  t. 


Then 


[-  kn  +  bt]  -  1,  if  t  <  t. 


M(t) 

Px(t  >  t>  -  ^  gx(j;  t). 

j-m(t)+l 


Let  Fr  designate  the  CDF  of  r^,  then 


F(_X)  (t)  -  1  -  F(X)(t)  -  Px{r  >  t>,  0<t. 


Furthermore,  from  the  Markovian  property  of  the  Poisson  process,  for  each 
(2.10)  F™  («)  -  f'rM 

M<Vl> 

+  2Z  8X(*  ti-l)  l1  -  Posf^tj.j)  +  1  -  j;  Mt  -  tj.j))] 
j-mft^jj+l 

This  formula  holds  for  i  =  1  too,  if  we  define  gx(j;  tg)  =>=  l{j  =  oj.  Notice  that 
(X) 

(t)  is  a  differentiable  function  of  t,  for  t^j  <  t  <  tj  (i  *=  1.  2,...,  ks  -r  kg). 
(X) 

Thus,  Fr  (t)  is  absolutely  continuous.  In  Section  4  we  develop  explicit  formulas 
(X)  (X) 

for  F^  (t)  and  FR  (t),  appropriate  for  Case  I. 


I 


Gv 

»A 

A 


£ 

Is 


ry 

A 

r 

!?: 

.2* 

I 

1 


I 
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3.  Linear  Boundaries  Crossing  Probabilities 

Let  h^  (i,  s)  designate  the  probability  that  a  Poisson  process,  with  intensity  X, 
assumes  the  value  i  at  time  s,  and  does  not  cross  the  linear  boundary  x(u)  =  bu  for 
0  <  u  £  s,  i.  e., 

(3.1)  h^(i;  s)  *  P^|sup(X(u)  —  bu  <  0;  0  <  u  <  s),  X(s)  =  ij. 

Let  1I~  be  the  region  below  x(u)  over  (0,  s],  and  tljjl'  the  region  above  it  (including 
the  boundary),  i.  e., 

(3.2)  «”«  {(u,  xY.  0  £  x  <  bu,  0  <  u  <;  s} 

and 

(3.3)  Ku»  *);  bu  £  x,  0  <  u  £  s). 


Conditioning  on  the  level  at  which  the  last  entrance  of  a  sample  path  to  <Ig  occurs, 
we  obtain  the  recursive  equation 


p(0;  Xs), 


if  i  -  0 


(3.4)  h,(i;  s) 


{ 


if  i  >  0,  s  £  tj 


a 

p(i;  Xs)  —  ^  ^  p(£;  Xt£)  h^(i  —  £;  s  —  t£),  if  i  >  0,  s  >  tj 

£-1 


Let  ^  =  X/b  and  let  H^(i;  s)  =  h^(j;  s).  Thus,  for  each  j=0,  1,  2,... 

j-0 


and  i  =  j  -+■  1,  j  +  2,. ..  we  obtain  the  recursive  relation 


{p(0;  iu),  j  =  v 

J 

Pos(j;  in)  -  p(£;  yu£)  Hf(j  -  £;  t^),  j  ^  1 

£-1 


j  =  0 


An  explicit  solution  of  (3.5)  can  be  obtained  by  using  the  theory  of  formal  power- 
series  (see  P.  Henrici  (1974);  p.  17). 
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For  6  £  (—1,  1)  and  6  —  0,  1,  2,...  define  the  power-series: 

hm  <•>  -  Z  W  eJ. 

j-0 

oc 

<3*6)  (fi)  -=  Posvj,  /Lt(j  -f  ff)) 

j-0 

oo 

px,tf  (0)  “  X  **  *0  +  «> 

j-0 

where  h^(0;  0)  m  1  end  p(0,  0)  m  1.  The  power-series  P^oU))  has  an  inverse, 

OQ 

given  by  Q*  (6)  -  q(X)  9n,  where 
A  — <  n 


i 

j  —  p(l;  u ),  n  —  1 
(-1)"  D „U),  n  ^  2 


where  Dn(/i)  is  the  determinant  of  the  n  x  n  matrix 


pCl;  u) 


p(2;  2 u) 
p(l;  u) 


(3.8)  An(/i) 


p(n;  n(A) 

p(n  —  1;  (n  —  1)m) 
p(n  —  2;  (n  —  2 )u) 
p(n  —  3;  (n  —  3)/u) 


p( 1;  u ) 


From  (3.5)  we  obtain  the  equation 


hm  (W  -  FM  <e)  -  hm  (S)  <p*0  (9)  -  "• 


or,  equivalently, 


(3.10) 


hM  m  p*o  (9> 


F*S  »>• 
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Hence, 


(3.11) 


H*  (0)  -  F*  (6)  Q*(0) 

X 


and,  for  each  i  *  0*  lf... 


(3.12)  Hx(i;  ti+tf)  =  23  <l(nX)  Pos(i  “  n?  (i  ~  n  +  *>M>. 

n-0 


In  e  similar  fashion  we  obtain  that  for  every  j  —  0, — ,  i  —  1  (1  —  1,  2,. ..) 

j 


(3.13) 


hx(j;  tj)  —  23  —  n;  (i  —  n)/x). 


n-0 


(X) 


Consider  a  linear  boundary  Xg(u)  —  —  £  +  bu,  for  an  integer  value  £.  Let  f ^  ^  (i) 


be  the  probability  that  the  process  crosses  this  boundary  for  the  first  time  at 
(t£+i,  i),  i  —  0,  1,  2r...  These  boundary  crossing  probabilities  satisfy  the  following 
recursive  equation  for  i  —  0,  1,  ... 


i-1 


(3.14)  f  (i)  —  p(i;  (£  +  i),u)  —  ifi  >  l}  f  .(j)  p(i  —  j;  (i  —  j)  ^)- 

J»0 


Define  the  power-senes,  for  0  €  (—  1,  1), 


(3.15) 


*  »  ' '  (X)  i 

fX,£“,-ZrA,£(')(’ 


1-0 


The  recursive  equation  (3.14)  yields  the  equation 


(3.16) 


f*  (6)  =  P*  (0)  (P*  (0)]'1, 
x,t  M  1  *»o  ’ 

where  *  =  Q*(0).  Hence,  we  obtain  by  multiplyjng  the  power-series, 

n  jU  X 


(3.17)  f  (i)  -  S'  q(M  pU  -  n;  (£  +  i  -  n )u\  f or  i  -  0,  1, 

n-0  " 


-7. 


13 


a 


21 


■J 


ft 


ft 


I 


r. 

u 


1** 


a 


A 


<  «J 


£l 
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4.  The  Distributions  of  Stopping  Times  in  Case  I 


The  continuation  region  (between  b^(t)  and  bR(t))  is  partitioned  into  four 


regions: 


Rj  =  {(t,  x):  0  <  t  £  tts-k1*  0  £  x  < 

^2  =  {^»  x^:  'kg-kj  <  x  £  'ko*  0  ^  x  <  bR(t)j; 

R3  —  {(t,  x):  kQ  <  t  £  •  b^(t)  <  x  <  bR(t)};  and 

R4  -  {(t,  x>.  tkQ+ki  <  t  £  tkg+ko,  bA(t)  <  x  <  bR(t)}. 


4.1  Distributions  over  R, 


Over  Rj,  r  —  tr,  while  Px{t  —  r A)  —  0.  From  the  Markovian  properties  we 
obtain,  for  i  =  1,...,  k8  —  kj,  j  —  0,  1,...,  kj  +  i  —  1, 


(4.1)  gx(j;  tj)  -  p(j;  ui) 


—  l{j  ^  kj  +  l}  p(k j  +  t\  iu)  hx(j  -  kj  -  tj  -  t£). 

£-1 


Hence,  for  1  <.  1  <;  ks  —  k,, 


(4.2)  Px{tr  >  tj  =  Pos(kj  +  i  —  1;  i,u) 

kj+i-1  j-kj 

-  >{■  *  -}  Z  Z  p,ki 

j-kj+1  £-1 

Changing  the  order  of  summation  we  obtain 


£;  in)  hx(j  -  kj  -  £;  tj  -  t^ 


(4.3)  Px{tr  >  tj  =  Pos(kj  ->-i  —  l;  i u) 


•  A 

!{i  ^  ‘}  ^2  p(kl  +  C;  iu)  Hx(l  -  1  -  ('•  Vt}- 
£-1 
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Thus,  for  every  i  —  1,...,  ks  —  kj,  we  obtain  from  (4.3)  and  (3.12), 


(4.4)  F  (tj)  =  1  —  Pos(k1  +  i  —  1;  iv) 


+  '{*  *  2> 


4*  tkf  a, 

•  Pos(i  —  1  —  £  —  n;  (i  —  £  —  n)u\ 

t-1  n-6" 


4.2  The  Distributions  over  R- 


Over  R2  we  obtain  similar  results  to  those  over  Rj,  with  only  slight  changes. 
Here,  the  CDF  of  is 

.  0,  i  -  k8  -  kj  4-  I,-  -,  kQ  -  1 

(4.5)  Fa<V-  (X) 

A  lf\  .  (0),  i  -  k*. 


The  CDF  of  rp  is  given  by 


(4.6)  (tj)  -*  1  —  Pos(ks  —  1;  in) 


ks-krl 


+  l{ks  >  kl  +  l}  P(kl  +  £;  iU)  HX(ks 

£-1 


-  k,  -  £  -  1;  t.  -) 


for  i  *=  ks  —  kj  4-  1,...,  Xq.  According  to  (3.12), 


H/ks  -  kj  -  £  -  1;  t5_,) 


ks-kj-£-l 


^  ]  q  k  Pos(ks  —  kj  —  £  —  n  —  1;  (i  —  £  —  n)*J- 
n-0 


4.3  The  Distributions  over  R- 


The  probability  distribution  (PDF)  of  over  Rg  is  given  by 


Px{tA  “k0  +  i}  =  fAk  (0’  1  =  1'—  kl’ 


y. 


§ 


r»:* 


9 


e 

Jl 


I 

I 


[fi- 


M 
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(X) 

where  f  (i)  is  given  by  (3.17). 

A‘k0 

The  probabilities  gjjj;  for  i  £  i  £  kj,  ere  given  by 

(4.9)  gx(j;  tkQ  +  i)-  l{l  £  j  £  kj 

(  p(j;  M(ko  +  i))  -f;  (€)  p(j  -  «  0 

l  A,Jt0 

+  l{kj  +  1  £  j  £  kt  -  l}  p(j;  M(kQ  4-  »)) 


-  «-) ) 


-  l{kj  +  1  £  j  £  kg  -  1} 

f(.X>.  (0  pfj  -  ti  M(i  -  «)) 

A,JC0 

-  l{kj  +  1  £  j  £  kg  -  l} 
j-kl 

X  *(kl  +  ■:  H  hx(j  -kt-m  tk  J 
■-1  0 

where,  according  to  (3.13),  for  kj  +  1  £  j  «£  kt  —  1, 

(4'10)  hx(j  ~  kl  -  ■?  ^ko+i-m) 

j-kj-m 

-  ^2  qnX)  p(j  —  kj  -  »  -  n;  u(k0  -hi  —  to  —  n)). 
n-0 


-10- 
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From  (4.9)  we  obtain,  for  i  £  i  £  kj, 


(X)  (X)  ,  ^ 

«•“>  Fd  1  -  f.  ‘V  -  Px{T  >  *i} 


1  —  F  (tj)  —  Pos(ks  —  1;  M(k0  4-  0)  +  Pos(i;  mOc0  4-  0) 


+  52  f  V .  («)  [Po»(ks  -  1  -  £;  u(i  -  £))  -  Poirti  -  £;  M(i 

1=0  A’ko 


+  l{k8  >  kj  4-  1} 


kg-k.-l 


X  P(kl  +  £;  €^}  Hx(k®  ~  kj  -  £  -  1;  tk  +i_J. 

£-1  0 


4.4  The  Distributions  over  Ri 


We  assume  that  ks  >  kj  +  1.  If  ks  —  kj  +  1  then  P^{r  £  “  1* 


Applying  the  Markovian  property  of  the  Poisson  process,  we  obtain  for 


1  £  »  <:  k8  —  kj  -  1, 


«•“>  P^(TA  -  V‘l  +  *>-  C„  <kl  +  0 


l 

“  22  gxlkl  €;  'kg+kjl 


/^\  | 

p(i  —  t\u\)  —  l|]  >  £}  y  ]  f  (m)  p(i  —  £  —  m;  u( i  —  £  —  m)N 

m-0  ’  ' 


where  g^(kj  4-  £;  tk^+k  )  is  given  by  (4.9). 


I.« 


I 


In  a  similar  manner  we  establish  that,  for  every  1  £  i  <.  ks 


J 

U3)  B\(kl  +  'kg+kj+i)  “  !{*  +  1  ^  j  ^  ks  “  kl  ~  *}  2^  *k(ki  +  £i  'kn+kJ 

VI  £_J  01 


f  i-l-1 

|p(j  —  t\  ^i)  —  l{£  <  i}  ^  "  f  (m)  p(j  —  £  —  in;  u( i  —  £  —  m)] 
1  m-0  A,t 


ence. 


k8-krl 


.14) 


Px{T  >  lk0+kj+i}  “  K1  <  ka  -  ki  -  3}  2  *xlkl  +  *  ^k/j+k.+i) 

j-i+1  1 

-  l{i  <  ks  -  kj  -  l} 


{  2Z  gx(kl  +  £;  'kg+kj)  [Pos(ks  -  kj  -  £  -  1;  iAt)  -  PosG  -  £;  i/i)] 

kg-kj-1 

£  gx(ki  +  *  VkJ  Po8(k8  “  kl  ~  €  “  15  *) 

£-i+l  0  1 


i-2  i-l-m 


lii  > l)  2  2Z  *x(ki  +  *  tk  k )  f(x) 

m-0  £-1  0  1  A,t 


(m) 


J"Pos(ks  —  kj  —  £  —  m  —  1;  u(\  —  £  —  m))—  Posii  —  t  —  m;  u(\  —  £  —  m)fjj. 


5.  Numerical  Example 

It  is  a  simple  matter  to  compute  the  CDF  of  the  stopping  times  (variables) 

tr,  ta  and  r  according  to  the  formula  given  in  Section  4.  In  Table  5.1  we  present 

these  distributions  for  Case  I,  with  kn  «=  15,  k,  =  15,  kB  =  28  and  b  *=  20. 

(X)  (X)  U  1  & 

We  present  F  (t),  F  (t)  and  F^  (t)  at  the  values  of  t4  «=  i/b 
(i  *=  1,  2,...,  k8  +  kg  —  1)  for  Poisson  processes  with  X  ■=  10  and  X  =  20.  Notice 
that  the  values  of  X  are  specified  for  a  convenient  time  unit,  which  is  traditionally 
take  to  be  the  MTBF  under  the  alternative  hypothesis,  and  a  specified  number  of 
systems  on  test  (see  Zacks  (1987)).  The  values  of  Table  5.1  were  computed  on  a  PC, 
using  a  TURBO-PASCAL  program,  which  is  available  upon  request. 
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Table  5.1:  Cunulative  Distributions  of  Stopping  Tlaes 

For  Truncated  Wald  SPRT  (kQ-^-15,  kg-28,  B-20) 


l 


X  -  10 

X  -  20 

| 

t 

rR(t) 

FA<t> 

F(t) 

FR(t) 

FACt) 

F  (t) 

| 

0.05 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

> 

? 

0.10 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.15 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.20 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

j 

0.25 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

| 

0.30 

0.000000 

0.000000 

0.000000 

0.000002 

0.000000 

0.000002 

\ 

0.35 

0.000000 

0.000000 

0.000000 

0.000005 

0.000000 

0.000005 

:■ 

0.40 

0.000000 

0.000000 

0.000000 

0.000013 

0.000000 

0.000013 

S 

0.45 

0.000000 

0.000000 

0.000000 

0.000028 

0.000000 

0.000028 

5 

0.50 

0.000000 

0.000000 

0.000000 

0.000056 

0.000000 

0.000056 

D 

0.55 

0.000000 

0.000000 

0.000000 

0.000100 

0.000000 

0.000100 

0.60 

0.000000 

0.000000 

0.000000 

0.000165 

0.000000 

0.000165 

0.65 

0.000000 

0.000000 

0.000000 

0.000258 

0.000000 

0.000258 

0.70 

0.000000 

0.000000 

0.000000 

0.000661 

0.000000 

0.000661 

oj 

0.75 

0.000000 

0.000553 

0.000553 

0.001728 

0.000000 

0.001729 

1* 

0.80 

0.000000 

0.003069 

0.003069 

0.004111 

0.000002 

0.004113 

0.85 

0.000000 

0.009555 

0.009555 

0.008836 

0.000007 

0.008844 

i 

0.90 

0.000000 

0.022050 

0.022050 

0.017319 

0.000020 

0.017339 

* 

0.95 

0.000001 

0.042105 

0.042106 

0.031269 

0.000044 

0.031312 

1.00 

0.000002 

0.070480 

0.070482 

0.052481 

0.000085 

0.052566 

1.05 

0.000005 

0.107088 

0.107094 

0.082541 

0.000149 

0.082690 

U 

1.10 

0.000013 

0.151118 

0.151131 

0.122503 

0.000243 

0.122747 

1 

1.15 

0.000027 

0.201242 

0.201269 

0.172631 

0.000373 

0.173004 

1.20 

0.000056 

0.255844 

0.255900 

0.232258 

0.000545 

0.232803 

1.25 

0.000109 

0.313230 

0.313340 

0.299814 

0.000764 

0.300578 

ft 

1.30 

0.000204 

0.371778 

0.371983 

0.372996 

0.001035 

0.374031 

r, 

f , 

1.35 

0.000367 

0.430045 

0.430412 

0.449057 

0.001362 

0.450420 

• 

1.40 

0.000635 

0.486825 

0.487460 

0.525136 

0.001749 

0.526885 

• 

1.45 

0.001061 

0.541171 

0.542232 

0.598567 

0.002198 

0.600766 

1.50 

0.001716 

0.592389 

0.594105 

0.667131 

0.002712 

0.669843 

v 

1.55 

0.002692 

0.640016 

0.642707 

0.729195 

0.003291 

0.732486 

1.60 

0.004105 

0.683788 

0.687893 

0.783771 

0.003936 

0.787707 

1.65 

0.015958 

0.713544 

0.729502 

0.830488 

0.004514 

0.835003 

$ 

1.70 

0.029219 

0.737448 

0.766667 

0.869395 

0.005117 

0.874512 

i 

1.75 

0.042745 

0.758349 

0.801093 

0.901008 

0.005788 

0.906796 

7\ 

1.80 

0.056490 

0.777129 

0.833618 

0.926085 

0.006545 

0.932630 

1.85 

0.070690 

0.794002 

0.864692 

0.945507 

0.007393 

0.952899 

>! 

1.90 

0.085617 

0.809013 

0.894631 

0.960181 

0.008329 

0.968510 

1.95 

0.101469 

0.822211 

0.923680 

0.970978 

0.009350 

0.980328 

{ 

2.00 

0.118164 

0.833688 

0.951852 

0.978670 

0.010451 

0.989121 

) 

2.05 

0.134752 

0.843573 

0.978325 

0.983895 

0.011627 

0.995522 

? 

2.10 

0.147980 

0.852020 

1.000000 

0.987127 

0.012873 

1 . 000000 

\ 

N. 

v" 

C 
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